Electronic transport through a quantum dot network 
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The conductance through a finite quantum dot network is studied as a function of inter-dot cou- 
pling. As the coupling is reduced, the system undergoes a transition from the antidot regime to 
the tight binding limit, where Coulomb resonances with on average increasing charging energies are 
observed. Percolation models are used to describe the conduction in the open and closed regime and 
contributions from different blockaded regions can be identified. A strong negative average mag- 
netoresistance in the Coulomb blockade regime is in good quantitative agreement with theoretical 
predictions for magnetotunneling between individual quantum dots. 



I. INTRODUCTION 

Arrays of insulating islands in two-dimensional elec- 
tron systems in the ballistic regime are often referred 
to as antidot latticesi^ Features observed in the magne- 
toresistance include commensurability peaks correspond- 
ing to quasi pinned orbits around 1,4,9, ... antidos^ and 
Aharonov-Bohm type oscillations superimposed on the 
first commensurability peak.?*S. The search for an artifi- 
cial band structure at B=0T is still ongoing. 

Here we present measurements on a finite array with 
a very small lattice constant of 120 nm, where the global 
electron density can be varied continuously with a metal- 
lic top gate. This allows us to monitor the transition from 
an antidot to a quantum dot array, that takes place when 
the electron density is reduced and the constrictions be- 
tween neighboring antidots enter the tunneling regime. 
Our sample is special in the sense that the extremely 
small lattice constant raises the charging energies to well 
observable levels and ensures that the lattice enters the 
tunneling regime well before the leads go insulating. 

In the following we primarily focus on electronic trans- 
port in the quantum dot network regime, which can be 
compared to conduction through granular or disordered 
materials. Related systems include arrays of metallic 
nano-crystalsiSiSiifliii, layers of semiconductor quantum 
dots^^, porous silicon^, 3D arrays of semiconductor 
nanocrystalsiiiiiii and organic molecular crystals^. How- 
ever, unlike most experiments on macroscopic samples, 
we are able to tune the inter dot coupling continuously 
and resolve individual Coulomb resonances due to the 
mesoscopic dimensions of our system. The local and 
global properties of the network are investigated by mea- 
suring across different terminals and phase coherence is 
probed by applying a perpendicular magnetic field. 



lithography to define the nanostructure under study. 
This patterning method relies on an atomic force micro- 
scope with a conducting tip to locally oxidize the surface 
of a GaAs heterostructure and thereby locally depicting 
the underlying 2DES (see Refs Il8l and Il9l for details). 
In this way a square lattice of 20 x 20 insulating islands 
with a lattice constant of a=120 nm was fabricated and 
enclosed by an insulating cavity with openings in the 
corners that serve as current and voltage leads (see inset 
Fig.l). Later a TiAu top gate was evaporated over the 
entire structure using a shadow mask technique. 



III. MEASUREMENTS 

At high electron densities clear commensurability 
peaks around 1 and 4 antidots appear with superimposed 
ballistic conductance fluctuations (see Fig. QJ, compara- 




II. SAMPLE 

Starting with a high quality GaAs/AlGaAs het- 
erostructure hosting a two-dimensional electron sys- 
tem (2DES) 34 nm below the surface, we used AFM- 



FIG. 1: Magnetoresistance measured from A to C across di- 
agonal 1 at different top gate voltages at T=90mK, commen- 
surability peaks around 1 and 4 antidots are marked. Inset: 
AFM-micrograph of the antidot lattice and the enclosing cav- 
ity. Bright regions are oxidized and correspond to depletion 
in the underlying two-dimensional electron system. 
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FIG. 2: (a) Conductance as a function of top gate voltage 
across diagonal l(A-C) and diagonal 2(B-D) at T=90mK. 
Bright regions are oxidized and correspond to depletion in 
the underlying 2DES. (b) Schematic of the lattice illustrat- 
ing quantum dot formation and the bond percolation model, 
(c) Double logarithmic plot of the averaged conductance over 
Akf — kf — fc/ pc with kfp C = kf at P c . The slope of the linear 
fit yields a value of £ CT = 1.2. 
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FIG. 3: Current as a function of top gate and bias voltage. 
White regions represent absolute current values above 5 pA. 
The upper two graphs are overview plots for diagonals 1 and 
2, dashed lines are guides to the eye and correspond to a 
scaling exponent of about 3 with respect to &f- The lower two 
graphs are blowups showing well defined Coulomb diamonds. 
All measurements were taken at a He bath temperature of 
90 mK at B = T. 



ble to measurements taken on a similar sample with a 
larger lattice constant by Schuster et a\&. This indi- 
cates that despite the small period, we still have a very 
symmetric two-dimensional potential modulation. It is 
worth noting, that pronounced commensurability max- 
ima in the magnetoresistance even occur for values ex- 
ceeding the resistance quantum. 

By applying suitable voltages to the top gate electrode, 
the electron sheet density in the 2DES can be tuned from 
about 2 to 5.5 x 10 15 m~ 2 . The carrier density in the 
lattice is about 1.5 x 10 15 m~ 2 lower than in the unpat- 
terned 2DES as determined from the Shubnikov-de Haas 
effect. Since the resistance in the leads is smaller than 
h/2e 2 down to top gate voltages below -100 mV it can 
be neglected in all measurements discussed here. Figure 
2(a) shows the conductivity measured across diagonal 1 
from corner A to C and across diagonal 2 from corner B 
to D as a function of top gate voltage at T=90 mK. As 
the voltage is lowered, the electron sheet density and the 
conductivity decrease until the Coulomb blockade regime 
is reached. This transition takes place at a top gate volt- 
age of about —50.8 mV, marked by V c in Fig. 2(a). To- 
wards even lower voltages a series of on average decreas- 



ing Coulomb peaks is observed until conduction com- 
pletely ceases below about —90 mV. It is worth pointing 
out, that the conductivity across both diagonals is very 
similar over the entire range of top gate voltages studied 
aside from mesoscopic fluctuations of the conductance 
caused by interference and interaction. This indicates 
the high symmetry and homogeneity of our sample. We 
would like to stress, that the physics in the ballistic anti- 
dot regime at high electron densities is in marked contrast 
to the quantum dot network regime close to and below 
V c . In particular the linear transport characteristics of 
an antidot lattice and the classical ballistic trajectories 
responsible for the commensurability maxima give way 
to nonlincarities owing to Coulomb charging and magne- 
totunncling between individual localized states. During 
this transition the conductivity changes by several orders 
of magnitude. 

In order to gain more insight into the electronic prop- 
erties of the Coulomb blockade regime, we measured the 
current as a function of top gate and bias voltage. As 
can be seen in Fig. 3, blockade is lifted at sufficiently 
high bias voltages and clear 'Coulomb diamonds' are re- 
solved. In contrast to analogous measurements on single 
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quantum dots, overlapping diamonds as well as streches 
in top gate voltage without Coulomb blockade are ob- 
served. This indicates the formation of a network with 
blockaded regions connected in series and in parallel. The 
charging energy for the individual resonances can be de- 
termined from the bias voltage maxima V ma x at the tips 
of the Coulomb diamonds according to: 



Echarging — cVrnax — ~t~ A^v, (1) 

where C is the capacitance and An is the quantum me- 
chanical single-particle energy spacing for the N th state. 
Neglecting A at, which makes a contribution of about 10% 
to the total energy, and applying a plate capacitor model, 
the area of the blockaded regions can be determined us- 
ing: 

C = e Q e Ga As C , B (2) 
a 

where Aqb is the area of the Coulomb blockaded re- 
gion and d is the distance between 2DES and top gate. 
This leads to an average size of about 30 unit cells 
(or dots) at V tg — —60 mV and about 4 unit cells at 
Vtg = —80 mV, if the diameter of the insulating discs is 
set to 80 nm based on the oxide height profile. 



are expected to be insulating, this is a reasonable as- 
sumption. This suggests a percolation transition with a 
characteristic scaling behavior of the typesSi 

a oc (Afcp)^, Ak F = k F (V tg ) - k F (V e ), (3) 

where £ is a critical exponent corresponding to a spe- 
cific quantity. From a double logarithmic plot of the 
conductance averaged across both diagonals as a func- 
tion of Akp, we extract a conductivity scaling exponent 
Ccr = 1.2±0.2 in the open regime (Fig. 2(c)). This can be 
compared with the calculated value of 1.32 ± 0.02 (Ref. 
Elf) for bond percolation in a classical square random 
resistor network, if phase coherence is neglected. Theo- 
retical work describing a similar scenario has also been 
set forth by Meir— in a model for the metal insulator 
transition in two dimensions. 

In the Coulomb blockade regime, we follow the 'links, 
nodes, and blobs model' introduced by StanlySi and 
Coniglic 24 . In this picture, the spanning network is de- 
composed into multiply connected 'blob bonds' and 'dan- 
gling bonds' forming 'blobs', that are in turn linked by 
individual 'cutting' bonds. In our system the blobs cor- 
respond to clusters of strongly coupled quantum dots 
that constitute the blockaded regions, while the 'cutting 
bonds' act as tunneling links. As the electron density is 
reduced, more and more QPCs pinch off and the clusters, 



IV. PERCOLATION ANALYSIS 

In the following we apply percolation theory2£ to an- 
alyze these findings. This can be done by considering 
the conductance of the entire lattice as being dominated 
by the constrictions between neighboring insulating discs 
forming quantum point contacts (QPCs). The area en- 
closed by four insulating islands can then be viewed as a 
quantum dot or artificial atom with four terminals con- 
necting it to its nearest neighbor quantum dots (see Fig. 
2(b)). As the electron sheet density is reduced, the con- 
ductivity of the QPCs decreases until the last channel 
pinches off and the QPC goes insulating. Owing to small 
inhomogeneities inherent to the fabrication process and 
the presence of stray background charges, the closing of 
the QPCs will be a statistical process. In a classical pic- 
ture, the QPCs can then be viewed as either conducting 
or broken resistors on a square bond-percolation lattice. 
Since the conductance staircase of a QPC can be approxi- 
mated to be linear with respect to the Fermi wavenumbcr 
kp = y / 27rn(V t g), where n(V tg ) is the top gate voltage 
dependent electron sheet density in the lattice, this pa- 
rameter is a natural choice for further analysis. In addi- 
tion we identify the percolation threshold P c for the open 
regime with V c = —50.8 mV and assume the fraction of 
conducting bonds, or QPCs, P, to be proportional to 
Akp = kp(V tg ) — kp(V c ). The exact relation depends on 
the width distribution of the QPCs, but close to the per- 
colation threshold P c = 0.5r& where half of the bonds 



diagonal 1 diagonal 2 

0.0 I(pA) l.o 0.0 ICpA) 1.0 




B(T) B CD 

FIG. 4: Conductance as a function of top gate voltage and 
magnetic field across diagonals 1 and 2 at a He bath temper- 
ature of 90 mK. Boxes highlight similar features associated 
with the same blockaded regions, while ovals mark features 
unique to one diagonal. 
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or blobs, shrink, leading to higher charging energies. The 
envelope functions in Fig. 3 correspond to a size scaling 
exponent of about (cb = 3 ± 1 as a function of — Akp 3H 
which is comparable to theoretical calculations for the 
blob size scaling exponent £b = 2.06 — 2.16 fRef l24T) and 
the mean cluster size 7 = 43/18 rj 2.4i211 

In order to make this interpretation more consistent, it 
might be more appropriate to consider a two stage pro- 
cess. First a QPC goes from open (a > 2e 2 /h) to tunnel- 
ing (cr < 2e 2 /h), before becoming practically insulating 
(a ~ 0). The first stage describes the transition from the 
open to the Coulomb blockaded regime, while the second 
one induces the transition from the Coulomb blockaded 
to an insulating state. The spanning cluster above fcf(K) 
then consists of open QPCs, while the cutting bonds in 
Coulomb blockade are formed by QPCs in the tunneling 
regime. In either case the interpretation of the critical 
exponents should be viewed as tentative owing to the ex- 
perimental uncertainty concerning the fraction P of open 
(tunneling) QPCs. We also point out, that a precise the- 
oretical understanding of Coulomb blockade scaling that 
includes effects from clusters connected in series and in 
parallel is still outstanding. 

Individual clusters of quantum dots in the network 
can be monitored by measuring the shift of the Coulomb 
peaks as a function of magnetic field across both diago- 
nals (Fig. 4). Since the magnetic field dependent energy 
variation of a quantum state is related to the exact shape 
and symmetry of its wave function, this variation can 
be regarded as a fingerprint of a specific cluster. Sim- 
ilar features in Fig. 4 can be attributed to the same 
cluster being traversed by current flow across both diag- 
onals, whereas differing features originate from clusters 
predominantly probed by transport across one of the two 
diagonals. This demonstrates, that transport is not dom- 
inated by a single quantum dot or small cluster close to 
one of the leads. 



V. PHASE COHERENCE 

Phase coherence across the spanning network was in- 
vestigated by measuring the magnetoconductance as a 
function of top gate voltage (Fig. 5). The averaged mag- 
netoconductance in the blockaded regime shows an in- 
crease with magnetic field by a factor of approximately 
3, which is significantly higher than the value of 4/3 pre- 
dicted by AlhassidS^ and measured by Folk et ali2& for 
single quantum dots. This increase occurs on a magnetic 
field scale that corresponds to one flux quantum per unit 
cell and is consistent with a magnetotunnling effect pro- 
posed by Raikh and GlazmannSi between elliptical 'elec- 
tron lakes'^ They predict a low field magnetoresistance 
of the form: 
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FIG. 5: (a) Conductance as a function of top gate voltage and 
magnetic field in the Coulomb blockade regime at T = 90 mK. 
A clear positive magnetoconductance for individual Coulomb 
peaks is observed, (b) Averaged magnetoconductance be- 
tween -54 mV and -70 mV in steps of 0.02 mV. Dashed lines 
mark a flux quantum through the unit cell (290 mT) and a 
fit using equation 0] (c) Conductance as a function of top 
gate voltage and magnetic field in the open regime, (d) Av- 
eraged magnetoconductance between 30 and mV in steps of 
0.25 mV. Dashed lines mark a flux quantum through the unit 
cell (290 mT). 



where Bo depends on the details of the tunnel bar- 
rier, but is typically of the order of h/e ■ l/d±d2, d\ and 
di being the semiaxes of the two electron lakes. A best 
fit to the magnetoconductance averaged with respect to 
top gate voltage in the Coulomb blockade regime (dashed 
curve in Fig. 5(b)) yields a value of 2.2 T for Bq. This 
corresponds to an average radius of about 40 nm, which 
is compatible with a single quantum dot confined to a 
unit cell. Evidence for this effect has also been reported 
by Voiskovskii and Pudalov- 8 , however without resolv- 
ing the increase in conductivity of individual Coulomb 
peaks. Since Ref. |2j] assumes perfect coherence within 
an electron lake and from the characteristic field scale of 
a flux quantum through a unit cell, we conclude, that the 
average dot radius can be considered as a lower bound 
for the phase coherence length in the Coulomb block- 
ade regime. Remnants of this effect are still visible in 
the open regime (Fig. 5(d)), but weak localization is 
more prominent. From the dip around ±5 mT we ex- 
tract a phase coherence length of about 300 nm in the 
open regime. 

In conclusion we have presented measurements on a 
multiply connected multi-terminal quantum dot network 
with tunable inter-dot coupling. Our sample is at a meso- 
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scopic scale, where collaboratively Coulomb blockaded 
regions can be discriminated and related to macroscopic 
properties. For strong coupling close to the percolation 
threshold, a classical random resistor network model with 
superimposed quantum fluctuations can be applied until 
charge quantization becomes important in the tunnel- 
ing regime. For weak coupling, Coulomb blockade dom- 
inates resulting in the theoretically predicted^SLi^i and 
experimentally observed 32 33 insulating state for T^O 
with current onset above a bias voltage threshold and 
hopping transport at elevated temperatures. A strong 
parabolic decrease in average magnetoresistance in the 
Coulomb blockade regime around B=0 T, is in good 
quantitative agreement with theoretical predictions by 
Raikh and Glazman^i. These findings complement re- 
ports by Wiebe et al^, who performed scanning tunnel- 



ing experiments, Hani et ali^i, who obtained local poten- 
tial information, and by Eytan et al^ who did scanning 
near field optical microscopy on two-dimensional perco- 
lating systems. We are able to quantitatively distinguish 
two length scales that are of fundamental importance for 
transport in the weak coupling regime. The first rep- 
resents the decreasing area of the Coulomb blockaded 
clusters while the second one describes wavefunction lo- 
calization to within the individual unit cells. Our mea- 
surements can also serve as an intuitive picture for the 
formation of the so-called 'Coulomb gap'^i that opens 
up around the Fermi energy as a function of electron lo- 
calization and Coulomb interactions. 

We thank Y. Gefen, R. Leturcq, C. Marcus, and D. 
Pfannkuche for valuable discussions. This work was sup- 
ported by the Schweizerischer Nationalfonds. 
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